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NONLINEAR  INTEORAL  EQUATIONS  OP  RADIATIVE  TRANSFER 


1.  INTRODUCTION 

The  area  of  mathematical  physics  which  we  shall  discuss 
Is  somewhat  unique  compared  to  the  others  discussed  In  this 
seminar.  A  peculiar  feature  of  this  subject  Is  the  fact 
that  within  the  mathematical  model  of  certain  radiative 
transfer  problems  there  are  both  linear  and  nonlinear 
equations  which  determine  desired  functions.  As  opposed 
to  the  usual  situation,  the  linear  equations  are  exact 
and  not  merely  linear  approximations  to  the  nonlinear 
equations.  This  fortunate  situation  permits  us  to  obtain 
many  exact  results  about  the  nonlinear  equatlcxis  by  a 
study  of  the  linear  equations. 

We  shall  discuss  one  of  the  relatively  simple  problems 
In  the  theory  of  radiative  transfer;  that  of  determining 
quantitatively  the  Intensity  of  radiation  (neglecting 
polarization  and  redistribution  of  absorbed  radiation) 

In  a  homogeneous,  plane-parallel  atmosphere  that  Is  of 
finite  thickness  and  Infinite  lateral  extent,  and  that 
Is  exposed  to  monochromatic  radiation  Impinging  on  one  side 
of  the  atmosphere  from  a  given  direction.  This  also 
describes  the  transport  of  monoenergetlc  neutrons  In  a  sub- 
crltlcal  slab  reactor. 

The  mathematical  formulation  of  this  problem  as  a 
scalar  linear  Integro-dlfferentlal  equation  dates  back 
to  Schwarzchlld  [64,63]  and  has  been  extensively  studied 
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[Cf.  16,  21,  30,  33,  34,  36,  38,  39,  46,  4?,  78].  If  a 
steady  state  prevails,  boundary  conditions  on  the  face  of 
the  atmosphere  allow  the  conversion  to  a  linear  Integral 
equation.  Ihla  equation  is  solved  uniquely  by  a  series, 
the  n-th  term  of  which  gives  the  contribution  of  radiation 
which  has  been  scattered  n  times  [Cf.  34]. 

In  theory  then  this  problem  is  ccMUpletely  solved. 

*^0  actual  computation  of  quantities  by  means  of  the  series 
Is  practical,  however,  only  If  either  the  atmosphere  Is  thin 
or  absorption  Is  strong  for  each  scattering.  It  is  In  an 
attempt  to  find  more  practical  methods  of  conputatlon  that 
various  authors  have  formulated  nonlinear  equations. 

Most  astronomical  observations  give  only  radiation 
reflected  from  or  transmitted  through  an  atmosphere,  whereas 
the  above  mentioned  solution  to  the  transport  equation  gives 
a  solution  throughout  the  atmosphere.  This  led  Ambarzumlan 
first  to  the  formulation  of  nonllneeir  Integral  equations 
for  the  Intensity  of  radiation  reflected  from  a  semi— 

Infinite  atmosphere  [l],  These  equations  were  derived  by 
observing  that  the  problem  for  a  semi— Infinite  atmosphere 
is  Invariant  under  the  addition  of  a  slab  of  finite  thick¬ 
ness.  "Ihls  Idea  has  been  exploited  by  Ambarzumlan  [4] 
and  Sobolev  [68,  69],  by  Chandrasekhar  [22,  23]  under  the 
name  of  "Principles  of  Invariance, "  and  by  Prelsendorfer 
[58-60],  Bellman,  Kalaba,  and  Wing  [11—13],  Redheffer  [62,  63], 
Ueno  [75,  76],  and  others  under  the  title  of  "Invariant 
Imbedding."  Many  of  the  derivations  of  these  equations 
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have  been  heuristic  and  are  similar  to  the  standard  heur¬ 
istic  method  for  deriving  equations  in  mathematical  physics 
by  particle  counting,  being  novel  in  that  the  physical 
parameter  of  atmospheric  thickness  is  varied. 

It  is  easily  shown  that  these  nonlinear  equations  can 
be  derived  rigorously  from  the  transport  equation  [8,  9, 

16,  36,  49].  For  more  ccmiplicated  geometries  this  seems 
preferable  since  the  heuristic  particle  counting  has  lead 
to  errors  [9]. 

We  shall  sketch  below  the  derivation  of  some  of  these 
nonlinear  equations  from  the  transport  equation.  We  shall 
also  show  that  the  same  functions  can  be  obtained  as  solu¬ 
tions  to  linear  singular  integral  equations.  With  these 
linear  equations  it  is  possible  to  do  precise  existence 
and  uniqueness  studies. 

For  the  sake  of  simplicity  we  shall  restrict  ourselves 
to  the  case  of  isotropic  scattering.  In  this  case  some 
of  the  derivations  are  deceptively  simple.  The  results 
which  we  shall  give  can  be  extended  to  anisotropic  scattering, 
and  we  shall  give  references  where  appropriate.  For  certain 
well  known  results,  which  are  derived  in  a  new  manner,  the 
method  of  derivation  is  given  since  it  is  the  one  most 
easily  extended  to  anisotropic  scattering. 

We  do  not  attempt  to  give  any  part  of  the  theories  of 
polarization  [22],  of  inhomogeneous  media  [13,  19«  49,  74] 
or  of  other  geometries  (8,  13,  22,  59]  and  other  boundary 
value  problems  [23].  We  do  not  touch  on  the  problems  of 


neutron  reactor  theory,  which  fit  naturally  Into  this 
subject  matter  but  In  which  one  often  asks  different 
questions  than  those  considered  here.  Neither  do  we 
consider  various  numerical  methods  ootained  by  replacing 
the  linear  transport  equation  by  a  system  of  approximate 
equations  [4,  6,  7,  20,  22,  71], 

In  summary  then  our  objective  is  to  take  the  simplest 
problem  which  is  of  physical  significance  and  which  will 
serve  to  display  the  derivations  and  useful  role  of  certain 
fomulatlons  other  than  the  transport  equation.  Of  par¬ 
ticular  Importance  are  the  questions  of  existence, uniqueness, 
and  practicality  for  computational  purposes. 

We  give  a  list  of  references  which  is  extensive  but 
by  no  means  exhaustive.  Additional  references  can  be 
found  in  many  of  the  papers  and  books  cited,  and  no  slight 
is  meant  here  by  the  omission  of  references  to  any  particu¬ 


lar  author  or  work. 
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2.  THE  TRANSPORT  EQUATION 

We  shall  consider  a  homogeneous  plane-parallel 
atmosphere  of  finite  optical  thickness  t.  Optical  depth 
into  the  medium  is  measured  by  a  variable  x,  0  <  x  <  t, 
in  a  direction  normal  to  a  face  of  the  atmosphere.  We  let 
I(x,  li,  denote  intensity  of  radiation  at  depth  x  in 

a  direction  specified  by  9,  the  inclination  to  the 
positive  X  axis,  and  by  the  azimuth  referred  to  a 

fixed  axis  normal  to  the  x  axis.  Here  M  •  cos  6. 

The  notion  of  homogeneity  should  perhaps  be  clarified. 
An  actual  atmosphere  can  have  variations  in  the  density 
and  in  the  mass  absorption  and  scattering  coefficients  [22], 
Tbe  optical  depth  is  introduced  as  a  variable  to  remove 
these  variations  from  the  equation.  If  after  this  trans¬ 
formation  of  variables  the  law  of  local  scattering  is 
constant,  the  medium  is  said  to  be  homogeneous. 

The  interaction  of  radiation  with  the  medium  results 
in  absorption  as  well  as  the  scattering  of  some  of  this 
energy  into  other  directions.  We  shall  assume  that 
radiation  fields  do  not  interact  with  each  other,  that 
the  frequency  of  the  radiation  is  unchanged  upon  scattering, 
and  that  absorbed  radiation  is  not  reemltted. 

We  consider  a  monochromatic  parallel  beam  of  radiation 
of  unit  average  Intensity  per  unit  time  falling  on  the 
entire  face  of  the  atmosphere  at  x  =  0  from  some  direction 
0).  We  suppose  that  a  steady  state  exists.  Rather 
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than  formulate  the  incident  field  as  a  boundary  condition 
in  the  form  of  a  6-function  in  direction,  we  consider  the 
reduced  incident  field  as  a  source  throughout  the  medium. 

If  we  consider  isotropic  scattering,  the  transport  equation 
for  the  intensity  of  diffuse  radiation  is  independent  of 
<P  and  is  given  by  [22,  34], 

(2.1)  M  (x,  v)  +  I(x,  M,  1^)  ■  J(x,  V) , 

where  the  source  function  J  is  given  by 

1 

(2.2)  J(x,  v)  =  exp[-x/i^]  +  ^  /l(x,  a,  i^)da, 

^  -1 

0<oo<l,  -I  <  <  1  ,  0  <  V  <  1. 

If  we  consider  the  medium  surrounding  the  atmosphere  as  a 
vacuum,  then  no  diffuse  radiation  enters  the  faces  of  the 
atmosphere.  This  gives  the  boundary  conditions 

(2.3)  1(0,  v)  =  I(T,  -n,  I/)  =  0  ,  0  <  M  <  1. 

We  have  given  the  transport  equation  for  isotropic 
scattering.  The  results  to  be  given  below  have  extensions 
to  the  case  of  anisotropic  scattering  where  cjd  in  (2.2) 
is  replaced  under  the  integral  by  a  phase  function  which 
is  usually  assumed  to  be  given  by  a  finite  expansion  in 


Legendre  polynomials  of  the  cosine  of  the  angle  between 
incident  and  observational  directions  [c.f.  22,  34]. 

Using  the  boundary  conditions  (2.3),  we  have  from  (2.1) 

X 

(2.4)  I(x,  l^)  *  n  /  ®xpt  (y  ”  x)/il]j(y,  U)dy 

^  0 

for  0  <  M  <  1 ,  and 

T 

(2.5)  I(x,  i/)  •  -rni  /  (y  ~  x)/4]j(y,  i')dy 

X 

for  -1  <  |X  <  0. 

If  these  relations  are  inserted  in  (2.2)  we  obtain 
the  familiar  linear  integral  equation  [22,  34] 

1  T 

(2.6)  J(x,t')  «  exp[-x/i/]  +  ^  /  /  exp[-lx-yl/M]j(y,i^)dy^ 

0  0 

The  integral  operator  without  cn  is  known  as  the  truncated 
Hopf  operator  A  and  it  has  been  intensively  studied 
[20,  30,  34,  38,  39,  50,  78].  This  operator  is  usually 
written  in  terms  of  the  exponential  integral  Ej^(x),  but 
for  later  computational  purposes  we  prefer  the  above  form. 

For  the  physical  application  it  is  the  solution  of 

(2.6)  for  V  in  the  interval  0  ^  ^  1  that  is  of 

importance.  For  the  mathematical  analysis  that  follows 
it  is  important  to  consider  solutions  for  all  nonzero 
complex  values  of  v.  We  shall  state  the  following 


Legendre  polynomials  of  the  cosine  of  the  angle  between 
Incident  and  observational  directions  [c.f.  22,  34]. 

Using  the  boundary  conditions  (2.3),  we  have  from  (2.1) 

X 

(2.4)  I(x,  u)  m  ^  /  exp[  (y  -  x)/M]j(y,  t')dy 

^  0 

for  0  <  M  <  1 ,  and 

T 

(2.5)  I(x,  M,  v)  m  /  exp[  (y  -  x)/M]j(y,  v)dy 

for  -1  <  <  0. 

If  these  relations  are  Inserted  In  (2.2)  we  obtain 
the  familiar  linear  Integral  equation  [22,  34] 

1  T 

(2.6)  J(x,i^)  •  exp[-x/v]  ^  exp[-lx-yl/M] J(y,«^)dy^ 

The  Integral  operator  without  O)  Is  known  as  the  truncated 
Hopf  operator  A  and  It  has  been  Intensively  studied 
[20,  30,  34,  38,  39,  50,  78].  This  operator  Is  usually 
written  In  terms  of  the  exponential  Integral  (x) ,  but 
for  later  computational  purposes  we  prefer  the  above  form. 

For  the  physical  application  It  Is  the  solution  of 

(2.6)  for  V  In  the  Interval  0  <  v  <  1  that  Is  of 
Importance.  For  the  mathematical  analysis  that  follows 
It  Is  Important  to  consider  solutions  for  all  nonzero 
complex  values  of  v.  We  shall  state  the  following 


well-known  theorem  [l6,  34,  36]  which  has  an  appropriate 
extension  to  anisotropic  scattering. 


Theorem  1 .  For  0  <  oj  <  1  and  0  <  T  <  oo  the 
equation 

(2.7)  J  -  CD  A(J)  =  B 

has  a  unique  solution  for  any  bounded  function  B  given 

00 

(2.8)  J  =  Soo"a”(B)  . 

n=0 

If  for  I  I  >  0,  we  take 

(2.9)  B  =  exp[-x/t'], 

then  J  is  analytic  in  the  extended  domain  I  I  >  0 

and  real  valued  and  nonnegative  for  real  i/  ^  0 . 

This  is  the  basic  existence  and  uniqueness  result 
which  expresses  the  solution  in  familiar  terms  as  the  sum 
of  contributions  from  n-th  order  scatterings.  The  series 
(2.8)  converges  geometrically  with  a  ratio  equal  to  cd 
times  the  maximum  eigenvalue  A(t)  of  A  .  Many  estimates 
of  X  as  a  function  of  r  have  been  made  [c.f.  20,  50] . 
Since  X  is  asymptotic  to  1  as  r  goes  to  Infinity, 
e.g.,  X(4)  =  .9,  the  series  computation  is  practical 


only  for  small  T  or  for  small  cu. 

Our  objective  below  is  to  discuss  other  formulations 
that  replace  (2.7)  by  other  equations  and  replace  (2.8) 
by  a  computation  which  is  practical  for  all  0  <  T  <  oo 
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3.  NONLINEAR  EQUATIONS 

We  now  derive  certain  nonlinear  equations  for  deter¬ 
mining  I  and  J  on  the  faces  of  the  atmosphere.  One 
motivation  for  this  Is  that  the  computation  of  I 
throughout  the  atmosphere  represents  unnecessary  work  for 
many  applications  since  It  Is  only  the  reflected  or 
transmitted  Intensities  that  are  observed  (p.  555  of  [s]). 
Another  motivation  Is  the  fact  that  the  outward  intensities 
on  the  faces  of  the  atmosphere  represent  missing  boundary 
conditions  for  the  transport  equation.  If  these  boundary 
conditions  can  be  determined  then  it  is  a  simple  matter 
to  compute  intensities  interior. 

Many  of  the  derivations  of  nonlinear  equations  are 
heuristic  arguments  based  on  ’’Principles  of  Invariance” 

[3,  22]  or  the  ’’Principle  of  Invariant  Imbedding”  [l3,  76). 
We  shall  not  state  any  of  these  principles  but  instead 
sketch  a  derivation  of  the  nonlinear  equations  from  the 
transport  equation. 

We  first  define  certain  standard  functions.  The 
scattering  function  S  and  the  transmission  function  T 
are  defined  by  [22] 

S(M,  P)  -  MI(0,  -M,  i'), 

(3.1) 

T(M,  P)  =  MI(T,  tM,  V)  ,  0  <  M  <  1. 
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By  (2.4)  and  (2.5)  these  are  merely  Laplace  transforms  of 
the  function  J  which  Is  defined  to  be  zero  outside  the 
Interval  0  ^  ^  that  Is 

T 

S(M,  ~  ^  exp[-y/M] J(y,  i^)dy, 

(3.2)  ^ 

T(M,  1^)  =  ^  exp[-(T-y)/M]j(y,  l/)dy. 

These  functions  are  well  defined  In  the  complex 
M-plane,  I  M  I  >  0,  as  entire  functions  of  the  variable 
/i~^.  In  this  setting  they  are  dependent  since 

(3.3)  S(M,  V)  =  exp[-T/M)T(-4,  l^)  . 

In  view  of  our  previous  comment  on  analyticity  of  J  In 
V,  S  and  T  are  analytic  In  M  and  v  for  all 
I  M  I  >  0  and  1  i/  1  >  0. 

Following  Busbrldge  [l6],  we  define  X  and  Y 
functions  by 

X(M)  ■  J(0,  M), 

(3.4) 

Y(M)  =  J(T,  M). 

From  (2.6)  It  Is  obvious  that 


(3.5) 


J(T  -  X,  M)  =  expI-T/M] J(x,  -M) 
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so  that 

(3.6)  Y(M)  =  exp[-T/M]X(-M)  . 

Using  (3.2)  and  (3.4)  In  (2.6)  we  see  also  that 

1 

X(/i)  -  1  +  ^  ^  S(a,  M)  , 

(3.7)  ^ 

Y(M)  =  exp[-T/M]  +  5  ^  T(a,  M)  ^  . 

By  the  definition  of  S  and  T  and  the  self- 
adjolntness  of  the  operator  A  ,  we  have 

S(M,  u)  =  S(»/,  M), 

(3.8) 

T(/i,  u)  =  T(l^,  M). 

This  fact  is  also  true  for  anisotropic  scattering  [3,  22,  52], 
but  the  proof  Is  more  complicated  since  the  operator 
replacing  A  in  general  Is  not  self-adJolnt  [52]. 

A  differentiation  of  equation  (2.6)  gives  by  Theorem  1 
the  following  theorem  [l6,  22]. 

*nieorem  2.  Let  J  satisfy  (2.6).  Then  for 
0  <  X  <  T  we  have 
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(3.9)  (x,  M)  +  M‘^J(x,  M)  =  /  J(x,  a)  ^ 

/  J(T  -  X,  a)  |2  . 

i!  0  « 

ft  T 

For  0  <  X  <  T,  the  onesided  derivative  is  given  by 

1 

(3.10)  (x,  M)  «  I  Y(M)  ^  J(T  -  X,  a)  ^  . 

This  theorem  has  an  appropriate  extension  to 
anisotropic  scattering  [52]. 

If  we  multiply  (3.9)  by  exp(-x/i')  and  exp[ -  (r-x)/!^] 
and  integrate,  we  obtain  the  familiar  equations  [l6,  22] 

(^  +  p)S(/i,  u)  =  X(M)X(i/)  -  Y(M)Y(l/), 

(3.11) 

(^  -  ^)T(M,  i')  =  X(#x)Y(i/)  -  YWXiv). 

If  these  are  inserted  in  (3.7)  we  obtain  Chandrasekhar's 
equations  for  0  <  M  <  1, 

x(M)  =  1  +  a. , 

YW  =  exp[-r/4l  .  »  /  Y(M)X(p)  -  X(M)Y(i>)‘ 

^  0  fj.  -  V 


(3.12) 
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If  (3.10)  Is  substituted  in  (3.9)  and  the  Laplace 
transforms  taken,  we  obtain  the  familiar  equations  [22] 

+  •^)S(/J.,i/)  + 

1  1 

«  [1  +  ^  ^  S(M,<r)  ]  [1  +  I  ^  S(v,  ff)  ], 

(3.13) 

^  T(li,v)  +  II  iti,v) 

=  [1  >  ^  ^SOi.o)  lie-"/*'  .  ^  ^Tiu.a)  f  ]. 

Equation  (3.10)  is  also  useful  in  obtaining  an 
instability  result  discussed  in  Sec.  6.  Setting  x  -  0 
we  obtain  [22] 

1 

(3.14)  II  (M,  r)  =  Yin,  ’■)  I  ^  Y(a,  T)  ,  T  >  0. 

Using  this  and  (3.6)  we  get  [22] 

1 

(3.15)  II  (M,T)  +  i  Y(M,T)  =  X(M,T)  ^  ^  Y(a,T)  ,  T  >  0. 

Analogous  definitions  and  results  have  been  given  for 
anisotropic  scattering  [22,  52].  For  isotropic  scattering 
problems  in  other  geometries,  Bailey  [S]  and  Bailey  and 
Wing  [9]  have  given  similar  rigorous  derivations  of 
invariance  principles  from  the  transport  equation.  One 
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Interestlng  feature  of  this  work  is  the  disclosure  of  an 
error  in  previous  heuristic  derivations  [9l. 

Equations  (3.12)  and  (3.13)  are  the  nonlinear 
equations  which  have  been  most  used  for  numerical  calcula¬ 
tions  [14,  24,  25,  26,  27,  28,  41 ] .  They  represent 
necessary  conditions  only  on  the  solution  J  to  (2.6) 
[22].  Before  considering  the  important  questions  of 
uniqueness  and  stability,  we  turn  in  the  next  section  to 
the  derivation  of  certain  linear  equations  which  are  also 
necessary  conditions  on  S,  T,  X,  and  Y. 


4.  LINEAR  SIMQULAR  EQUATIONS 


We  want  now  to  derive  certain  linear  equations  as 
necessary  conditions  on  the  unique  solution  J  to  equation 
(2.6).  For  the  Isotropic  case  a  derivation  of  these 
equations  from  the  nonlinear  X  and  Y  equations  Is  due  to 
Busbrldge  [l6j  36].  (See  also  Crum  [29].)  Busbrldge's 
approach  has  never  been  extended  to  similar  equations  for 
anisotropic  scattering.  For  this  reason  we  give  a 
different  derivation  for  the  Isotropic  case  which  does 
readily  extend  to  the  anisotropic  case  [54]. 

We  apply  the  operator  A  to  equation  (2.6).  With  A 
expressed  as  in  (2.6)>  it  Is  an  easy  computation  to  show 
that  for  V  replaced  by  any  complex  number  z  outside  the 
Interval  [-1,  l]  we  have  the  equation  for  A  (J) 

1 

(4.1)  uuA  (J)(x,z)  -  (a»A)^(j)(x,z)  =  ^  dxp[-x/z]  J 

-1 

1 

-  zw  r  exp[-x/a)  ^ 

2  J  z-c 
o 

-  expt-xA)  J  do  . 


By  the  uniqueness  result  of  Theorem  1,  we  obtain  the 
solution 
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(4.2)  (dA  (J)(x,z) 


1 


J{X,^)  ^  J 

-1 


do  _  sou 
z-o  ”  T 


J 


z-a 


o 


do 


-  exp[-T/z] 


The  last  term  follows  from  (3.5). 

Combining  this  with  ^2.6)  we  obtain  linear  conditions 

on  J. 

Theorem  3.  For  complex  z  outside  the  Interval  [-1,  l] 
the  unique  solution  to  (2.6)  also  satisfies  the  equation 

1 

(4.3)  X(z)j(x,z)  =  exp[-x/z]  -  ^  f  do 

o 

-  exp(-TA)  J  do  . 

o 

The  function  X,  Is  defined  by 

(4.«)  X(z).1.2f/^, 

-1 

and  X  la  merely  a  parameter  satisfying  0  <  x  <  t  . 

If  J(xp)i)  and  J(t-x,u)  are  known  for  0  <  p  <  1/ 
then  (4.3)  defines  their  meromorphlc  extension  to 
complex  z  outside  [-1,  1].  If  J  Is  to  satisfy  (2.6),  It 
must  be  analytic  In  all  z  for  tz|  >0.  This  means  that 
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the  right  hand  aide  of  (4.3)  must  vanish  at  the  zeros 
of  X. 

The  function  X  Is  even  and  has  only  two  real  roots. 
For  0  <  tt>  <  1  these  are  given  by  -  1/k  where  k  Is  the 
nonzero  root  of 

(4.5)  2  k  -  (B  to  . 

For  0)  «  1,  X  has  a  double  root  at  Infinity. 

We  demand  that  the  right  hand  side  of  (4.3)  vanish 
to  the  same  order  as  X.  For  0  <  w  <  1  we  obtain  the 
conditions 


(4.6)  exp[±  kx] 


J 


o 


1  i  ko 


(V 


exp[t  kr) 


I 


o 


J(T-X,g) 

1  +  kc 


For  IB  -  1,  we  have 
1 

2  -  J  [j(x,g)  +  J(t-x,o)J  do  , 

0 

(4.7) 

1  1 

T  J  j(T-x,o)do  -  2x  +  J  o[j(x,o)  -  J(T-X,o)]d0 
o  o 

We  have  therefore  obtained  two  linear  constraints  on  J 
for  each  x  by  demanding  analytlclty  In  the  variable  z. 


do  . 
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Now  we  can  use  PlemelJ’s  formula  [57]  to  obtain 
singular  equations  for  determining  J(x,u)  and  J(t-x,u) 
on  0  <  M  <  1.  For  u  In  0  <  m  <  1  we  take  limits  of  (4*3) 

In  the  complex  plane  as  2  tends  to  u  first  with  Im(z)  >  0 
and  then  with  Im(z)  <  0.  Adding  these  limits  and  using 
analytlclty  of  J,  we  obtain  the  singular  equation  [54] 

1 

(4.8)  Xq(u)J(x,u)  =  exp[-x/M]  -  J  do 

o 

.  M|  exp[-VM)  f  d,  . 

o 

The  singular  Integrals  are  computed  as  Cauchy  principle 
values.  The  function  is  given  by 

(4.9)  =  1  +  ^  -tn  . 

If  we  evaluate  (4.8)  at  x  =  0  and  at  x  =  t  and 
recall  (3*4)  we  obtain  the  equations  [16,  29j  31i  53] 

X^Mx(u)  ^  exp[-T/ul  /  2^ 

0  0 

V(p)Y(u)  -  exp[-Vu)[l  -  f  f  ^ 

o  o 


(4.10) 
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If  equation  (4.3)  le  miltlplled  by  exp[-x/v]  and 
by  eXp[-(T-.x)/v]  and  Integrated  on  x,  we  obtain  the 
following  analogous  equations  for  3  and  T  functions  [52]. 
The  linear  singular  equations  are 


(4.11) 


.*p[-r(i  .  i)D  -  1*5  J 


-  lig  .xp[-T/ul  J 


(expt-T/ul  -  exp[-T/vl^  "  ^  J 


-  iif  .xp[-T/u)  J  SK|1S£ 


For  0  <  (I)  <  1  the  linear  constraints  are 


(4.12)  (i  _  exp[-  T(i  i  k)]) 

1  -  KV 


S(vja)da 
1  +  ko 


+  ■J  exp(+  kr] 


J 


o 


T(vjO)da 

1  i  ko 


For  0)  ■■  1,  the  linear  constraints  are 


-22- 


1 

2v(l  -  exp[-T/v])  »  J  [s(v,a)  +  T(v,o)J  da 

o 

1 

(4.13)  T  J  T(v,o)do  »  2v^  -  (1  +  t/v)  exp[-T/v]^ 

0 

1 

+  J  o[s(v,o)  -  T(v,a)]  da 
o 

In  view  of  (3.11)  we  expect  solutions  to  (4.11)  to  be 
given  simply  by  solutions  to  (4.10). 

Even  though  (4.10)  is  Important  in  l>;s  relation  to 
scattering  and  transmission  functions,  (4.8)  for 
X  and  T-x  are  equally  Important  since  they  determine  the 
source  function  J  at  any  symmetric  depths  in  the  slab. 

Any  numerical  method  for  treating  (4.10)  applies  equally 
well  to  (4.8).  This  is  of  particular  significance  In 
reactor  criticality  problems  [3?].  In  a  paper  currently 
being  written  with  Mr.  Leonard,  we  are  showing  that  a 
study  of  (4.8)  leads  both  to  the  criticality  condition 
for  slabs  and  spheres,  l.e.,  the  functional  relation 
between  t  and  the  average  neutron  production  per  fission, 
and  to  Fredholm  equations  for  the  neutron  density  In  these 
critical  reactors.  This  is  related  to,  but  simpler  than, 
some  current  work  In  neutron  transport  [21,  45,  48,  8o]. 

Although  (4.3)  and  the  resulting  equations  have  been 
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derived  from  the  transport  equatlon«  they  express  some 
physical  principle  of  duality.  The  transport  equation, 
which  expresses  differential  changes  In  I  Interior  to  the 
slab,  has  been  replaced  by  an  equation  which  relates  J, 
and  hence  also  I,  at  x  and  t-x  by  means  of  the  functional 
dependence  of  J  on  variations  of  the  angle  of  Incident 
radiation  on  the  face  of  the  ataK}sphere.  Ve  refrain  from 
designating  this  as  a  new  principle. 
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5.  THE  UNIQUENESS  PROBLEM 

Before  considering  the  merits  of  various  equations  for 
numerical  computations  it  is  inoperative  to  determine  whether 
or  not  each  set  of  equations  has  a  unique  solution.  For 
some  reason  this  problem  is  often  confused  with  the 
problem  of  xiniqueness  in  the  physical  problem  supposedly 
represented  by  the  mathematical  model. 

We  Illustrate  with  a  quote  from  Ambarzumian  [3,  p.  539], 
where  he  is  discussing  the  symmetry  of  the  function  S 
(p  in  his  notation)  for  a  semi-infinite  atmosphere 
"...  However,  since  equation  (33.24)  can  have,  from  its 
physical  significance,  only  one  solution,  the  function  p(ri,0 
must  be  symmetrical:  ...". 

The  conclusion  about  symmetry  reached  by  this  argument 
is  correct.  However,  the  reasoning  is  incorrect  since, 
except  in  the  conservative  case  of  lu  *  1,  the  equation 
he  is  discussing  has  two  symmetric  solutions  expressed  in 
terms  of  the  two  solutions  [I6]  to  Ambarzumian 's  equation 
(eqn.  (33»28)  on  p.  559  of  [3I)  or,  what  is  the  same  thing, 
ChandraselChar's  H  equation. 

We  quote  from  Chandrasekhar  [22,  p.  208]  as  an 
illustration  of  a  more  appropriate  viewpoint.  "The 
reference  to  the  question  of  uniqueness  of  the  solution 
in  the  preceding  paragraph  draws  attention  to  a  fact 
which  we  have  so  far  Ignored,  namely,  that  there  are. 

Indeed,  no  grounds  for  believing  that  the  principles  of 
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Invarlance,  by  themselves^  can  lead  to  determine  solutions 
of  the  various  problems.  The  invariances  pnd  the  equations 
they  lead  to  are  only  necessary  conditions;  it  Is  by  no 
means  obvious  that  they  are  also  sufficient . ” 

The  confusion  of  physical  models  with  mathematical 
models  often  causes  the  mathematician's  query  to  be 
termed  "academic”  when  he  asks,  "But  do  the  equations  have 
a  unique  solution?"  With  regard  to  some  of  the  equations 
given  In  previous  sections  we  shall  see  Just  how  relevant 
Is  this  question. 

Let  us  consider  first  the  nonlinear  equations  (3.12) 
for  X  and  Y  functions.  Chandrasekhar  realized  the 
Importance  of  the  uniqueness  problem  for  these  equations, 
and.  In  fact,  he  shows  [22]  that  these  equations  have  a 
one  parameter  family  of  solutions  for  uu  «  1.  The  unique¬ 
ness  problem  Is  also  considered  by  Busbrldge  [16].  We.[5l] 

recently  gave  a  complete  discussion  of  the  solutions  to 
equations  similar  to  (3*12)  but  with  ou  replaced  under  the 
Integral  by  a  nonnegative  "characteristic  function"  whose 
Integral  Is  bounded  by  1/2.  We  have  also  studied  the 
uniqueness  question  for  similar  equations  In  anisotropic 
scattering  [52]. 

For  isotropic  scattering  we  show  In  [51I  that  for 
u)  =  1  Chandrasekhar  found  only  half  of  the  solutions 
to  (3.12).  We  designate  by  and  the  functions  (3*^)* 
Then  all  solutions  to  (3«12)  for  'X'  -  1  are  given  by 
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X(u)  -  (1  +  bYu)X^(M)  +  |i(a+bu)[xQ(u)  +  Yq(m)]  , 

(5.1) 


Y(u)  -  (1  -  bYM)Yo(u)  -  u(a  -  bu)[xQ(M)  +  Y^{u)]  . 


The  constant  y  Is  given  by 


(5.2) 


Y 


^1  yi 


with 


1 

(5.3)  *1  -  5  J  Xq(v)v^  dv  , 

o 


The  parameters  a  and  b  are  arbitrary  solutions  to 


(5.4)  b(bY^  +  2a  y-2)  -  0. 

The  solutions  with  b  -  0  are  given  by  Chandrasekhar  [22]. 
For  0  <  «B  <  1  the  solutions  to  (3.12)  are  given 


by 
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X(li)  .  fl  +  (f-gP)l^lIM^(f»*-BP)  ]  x„(u) 

L  1  -  (kur 


1  -  (iC4) 


(5.5) 


Y(u)  .  [l  -  i/'VrgP)n-|^^^fyBP)]  y„(h) 


.  (fS-ga)u-ku^(fBtgtt)  y.  ,  . 
1  -  ° 


The  constant  k  Is  given  hy  (4.3)  and  constants  a  and  p  are 
given  by 


-i-fj 


*„(v) 


"5  J  1+kv 
o 


dv  , 


(5.6) 


p-fl 


^  *„(v) 

1+kv 


dv  . 


The  parameters  f  and  g  are  arbitrary  points  on  the 
hyperbola 


(f^-g^)(a^-P^)  -  2k(fofgp).. 


(5.7) 
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These  results  were  easily  obtained  by  first 
describing  all  solutions  to  the  linear  equations  (4.10). 
The  totality  of  such  solutions  Is  given  either  by  (5.1) 
and  with  arbitrary  parameters  a  and  b  or  by  (3.5)  and 
arbitrary  parameters  f  and  g.  The  constraints  (3.3)  on  a 
and  b  and  (5.6)  on  f  and  g  are  the  only  additional 
restriction  Imposed  by  the  nonlinear  equations  (3. 12). 

So  we  have  a  set  of  equations  (3*12)  which  were  first 
derived  from  heuristic  "Invariance  Principles".  If  the 
variable  u  Is  restricted  to  the  values  0  <  u  <  1  of 
physical  lnterest«  these  equations  alone  do  not  specify 
a  unique  solution.  It  is  easily  shown  that  the  linear 
constraints  (4.6),  evaluated  at  x  >  0,  are  additional 
equations  needed  to  select  the  desired  functionals, 

X  and  Y,  of  the  solution  J  to  (2.6).  But  it  Is  hardly 
likely  that  these  constraints  would  be  derived  by  a 
heuristic  particle  counting.  They  come  from  the 
mathematical  condition  of  analytlclty  In  a  variable 
which  Is  usually  restricted  to  the  Interval  0  <  u  <  1 
as  the  cosine  of  an  angle  in  the  Interval  [O, 

The  fact  that  equations  (3«12)  do  not  have  a  unique 
solution  does  not  mean  that  they  are  wrong.  It  Just  means 
that  they  do  not  constitute  a  complete  set  of  equations. 
Without  a  mathematical  Investigation  of  uniqueness  of 
solution  to  a  set  of  equations  It  Is  not  always  possible 
to  say  that  the  some  needed  equations  may  not  have  been 
overlooked. 
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As  we  have  already  remarked,  the  linear  equations 
(4.10)  have  a  two  parameter  family  of  solutions.  Ve  know, 
however,  that  the  constraints  (4.6)  are  also  necessary 
conditions.  These  constraints  and  equations  (4.10)  are  a 
complete  set  of  equations  with  a  unique  solution.  We 
shall  see  In  Section  7  that  these  equations  can  be 
transformed  to  a  form  Ideal  for  numerical  computation. 

We  have  not  discussed  the  uniqueness  problem  for  the 
nonlinear  Integro-dlfferentlal  equations  (3.13).  We  shall 
do  this  along  with  other  considerations  In  the  next  section. 
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6.  INSTABILITY  RESULTS 

In  this  section  we  wish  to  discuss  the  appropriateness 
of  the  various  nonlinear  equations  for  numerical  computations. 
We  are  particularly  Interested  In  the  problem  of  stability. 
Since  solutions  to  equations  are  usually  approximated  by 
nxmerlcal  solutions  to  approximate  equations,  it  Is  helpful. 
In  bounding  errors,  to  know  properties  of  solutions  to  the 
perturbed  equations.  In  particular,  their  degree  of 
approximation  to  the  solution  to  the  exact  equations. 

We  wish  first  of  all  to  consider  the  nonlinear  X 
and  Y  equations  (3«12).  To  understand  the  Importance 
for  numerical  computations  of  the  multiplicity  of  solutions 
given  In  Sec.  5  we  need  a  result  about  the  function  J. 

Theorem  4.  Let  e(T)  denote  the  maximum  eigenvalue 
of  the  operator  JL  .  Then  e(T)  <  1  for  0  <  t  <  •, 
and  the  solution  J  to  (2.6)  Is  an  analytic  function  of 
ou  for  all  complex  lu  In  the  domain  [uul  <  l/e(T). 

This  theorem  Is  easily  established  by  standard 
spectral  theory  arguments  about  the  inverse  (I  —  uu  JL  )  ^ 
for  the  compact,  positive  definite,  self-iadjoint  operator 
JL  . 

This  result  shows  that  the  functions  and 

have  power  series  expansions  In  uu  which  converge  In 

the  circle  |uu|<  l/e(T).  We  now  show  that  this  falls  to 

be  true  for  all  solutions  (5.1)  and  (5.5)  to  (3.12)  other 

than  X„  and  Y„ . 
o  o 


The  function  (5*5)  contain  the  conatant  k.  As  a 
function  of  uj«  k  Is  determined  by  (4.5)«  which  we 
rewrite  as 


(6.1) 


i^2n-2 


•••]»  |h|<  1. 


It  Is  easily  seen  that  k  ■  1  at  (u  -  0>  but  that  k 
as  a  function  of  complex  w  has  a  singularity  at  u)  -  0 
and  a  branch  point  of  order  2  at  u)  ■■  1.  Therefore  k 
In  (5.5)  Is  not  an  analytic  function  of  uu  even  In  the 
domain  |u)|  <1.  The  functions  given  In  (5.1)  are 
appropriate  limits  of  those  given  In  (5.5)  as  w  tends 
to  1  and  k  tends  to  0«  where  k  has  a  branch  point 
as  a  function  of  ou. 

We  obtain  then  another  characterization  of  the 

functions  X.  and  Y„. 

o  o 

Theorem  5.  There  Is  a  \inlque  solution  to  equations 

(3.12)  \hlch  Is  analytic  In  lu  for  |tul<  l/e(T),  namely 

the  functions  X_  and  Y.. 

— ■  ■  -  o  — —  o 

A  standard  numerical  technique  for  solving  functional 
equations  such  as  (3.12)  Is  by  Iteration  [3*  20,  26,  27, 
32,  61].  Let  the  Iterations  be  defined  by 
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*n+l(“) 

(6.2) 


,  ,  u»  r  . 

^  +  T"  J - \rr^ - 

0 

expi^/u, .  j 


This  generates  the  n-th  partial  sums  in  the  power  series 
expansion  of  the  solutions  and  This  procedure  is 

also  equivalent  to  introducing  in  (3.12)  the  power  series 
expansions  in  powers  of  m  and  equating  coefficients,  a 
standard  perturbation  analysis. 

By  Theorem  5  this  iterative  method,  if  computed  with 
absolute  precision,  leads  to  the  desired  solution  X 

o 

and  Yq  to  (3.12)  by  reason  of  the  analyticity.  But 
what  of  the  many  soluticms  (5.1)  or  (5.5)?  These  are 
solutions  to  (3.12)  \diich  are  not  analytic  in  lu  and  in 
theory  will  not  be  reached  by  the  iterative  method  as  long 
as  the  initial  functions  in  the  iteration  are  not  chosen 
to  be  one  of  these  extraneous  solutions. 

There  are,  however,  extraneous  solutions  to  (3.12) 
\diich  are  arbitrarily  close  to  the  desired  solution 
(Xo»  Yq)*  It  is  to  be  expected,  as  has  been  observed, 
that  as  the  iterations  seem  to  be  converging,  numerical 
errors  introduce  the  Influence  of  the  extraneous  solutions 
and  cause  oscillations.  So  in  theory  the  iterative  methods 
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wlll  converge  to  a  xmlque  solution,  but  in  practice  the 
presence  of  extraneous  solutions  arbitrarily  near  the 
limit  Introduces  Instability  In  the  calculations.  Only 
by  extreme  accuracy  can  this  Instability  be  avoided  In 
treating  the  X  and  Y  equations  (3.12). 

The  nonlinear  eqtiatlons  (3.12)  are  obtained  from 
(3.9)  without  any  use  of  (3.10).  This  results  In 
nonuniqueness  of  solutions.  Rather  than  consider  (3.9) 
and  (3.10)  we  now  consider  the  equation  (3.13)  for  the 
function  S.  This  Is  a  nonlinear  Integro-dlfferentlal 
eqxiatlon  for  determining  S  as  a  function  of  the  total 
optical  thickness  T  of  the  atmosphere.  The  Initial 
condition  on  S  for  t  >  0,  l.e.,  no  atmosphere.  Is 
naturally 

(6.3)  S(u,  |i^,  0)  ■  0  . 

Equations  (3.13)  and  (6.3)  represent  the  conversion 
of  a  boundary  value  problem  for  linear  equations  (2.1) 
and  (2.3)  to  an  Initial  value  problem  for  Rlcattl  type 
equations.  The  computational  advantage  of  such  a 
formulation  Is  obvious  and  has  been  much  emphasized 
[13,  63].  Computations  with  these  equations  have  been 
performed  for  Isotropic  scattering  [14].  . 

If  we  specify  an  Initial  function  for  S  at  t  - 
by 

(6.4)  S(^,  V,  Tq)  -  s(u,  v)  , 
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we  can  readily  convert  (3*13)  to  a  nonlinear  Integral 
equation  of  Volterra  type 

s(u,  V,  t)  -  ^^-JsCu^v)  +  J  exp[=ii^|^]x 

(6.5) 

1  1 

^1  +  y  J  S(|i>o,t)~J^l  +  y  J  S(v,0>t)^Jdt. 

0  0 

One  can  show  that  for  sufficiently  small  Bxtd  s 

the  right  hand  side  of  (6.5)  Is  a  contracting  operator 
and,  hence,  that  the  equation  has  a  unique  solution. 

This  solution  can  be  continued  for  larger  values  of  t. 

In  recent  reports  [55,  56]  we  have  studied  solutions 
to  (3.13)  and  a  family  of  Initial  conditions  (6.4).  This 
Is  actually  equivalent  to  studying  equations  (3.14)  and 
(3.15).  Of  particular  Interest  are  a  family  of  singular 
solutions,  l.e..  Infinite  for  finite  T,  and  their  nearness 
of  approach  to  the  desired  solution  for  Increasing  r. 

We  designate  the  solution  to  (3.13)  and  (6.3)  by 
Sq.  It  can  be  shown  that 

(6.6)  1^  Sq(u,v,t)  -  ^  H(n)H(v), 

where  H  Is  Chandrasekhar's  notation  (see  also  Ambarzumlan 
[3])  for  the  limiting  value  of  the  function  as  t 

tends  to  Infinity. 


» 
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We  consider  first  the  conservative  case,  l.e., 
u)  ■  1.  It  can  be  shovn  [51,  53,  55]  that  as  t  tends 
to  Infinity 

(6.7)  Sq(u,v,t)  -  H(u)H(v)  +  0(i). 

So  Sq  approaches  the  limit,  given  In  (6.6),  only  as 
1/t. 

By  use  of  (3.14)  and  (3.15)  and  (5.1)  —  (5.3)  It 
Is  possible  to  exhibit  solutions  to  (3.13)  and  Initial 
conditions 

(6.8)  S(^,v,0)  -  4(1  +  2  Buv)|iv  . 

We  have  given  explicit  solutions  for  (6. 8) [55],  but  we 
only  give  them  here  In  the  form 

S(h,v,t)  -  Sq(w,v,t)  +  0(i) 

B 

+ - - —  0(1) 

l-2Bj(x^(t)  +yi(t)2dt 
0 

For  any  B  >  0  these  solutions  are  Infinite  for  some 
finite  value  of  t.  By  choosing  B  sufficiently  small 
and  positive  and  t  sufficiently  large  It  Is  possible 
to  find  a  singular  solution  which  comes  arbitrarily  close 
to  the  solution  before  It  becomes  Infinite. 
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For  0  <  uj  <  1  there  are  also  8lng;ular  solutions 
to  (3»13)«  but  then  these  can  come  within  a  distance 
of  Sq  proportional  to  the  constant  k  and  no  closer. 
Since  the  Initial  conditions  and  solutions  are  more 
complicated  In  this  case«  we  merely  refer  to  [56]  for 
details. 

The  nonlinear  Integrodlfferentlal  equations  with 
Initial  conditions  are  attractive  from  the  numerical 
viewpoint.  The  above  results  show«  however^  that  there 
Is  a  limit  to  the  accuracy  which  can  be  achieved  for 
large  values  of  at  least  In  the  conservative  and 
near-conservative  cases. 
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7.  LINKAR  PRBDHOm  BQHJATIONS 

lito  conclude  trlth  a  quick  survey  of  some  results 
obtained  from  the  linear  singular  eqviations  (4.10)  and 
linear  constraints  (4.6)  or  (4.7).  nie  method  for 
treating  these  equations  is  to  Invert  the  singular 
operator.  By  analytic  continuation  the  transformed 
equations  can  then  be  transfoxmed  to  Fredholm  equations 
with  continuous  kernels,  nie  seme  program  has  been 
carried  out  for  anisotropic  scattering  [ 34] «  but  in  less 
complete  fom. 

For  siiig>llcity  we  merely  give  here  the  results  for 
m  >  1.  Conqplete  details  are  given  in  [53]  for  X  and 
Y  equations  with  a  characteristic  function  f  instead 
of  the  constant  a>«  and  in  [34]  for  anisotropic 
scattering.  The  equations  for  X  and  Y  functions 
have  been  programned. 

By  consistent  use  of  the  theory  of  singular  integral 
equations  and  analytic  function  theory,  we  have  obtained 
the  following  formulae. 

We  designate  by  H  the  limit  of  the  Xq  function 
as  T  tends  to  infinity.  Then  H  is  given  on  0  ^  M-  ^  1 
by  the  foxvula  [33] 

(7.1)  H(n)  -  (l-Hi)  J  dt] 

o 
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The  fmctlon  \q  is  given  by  (4.9)  for  m  ■  1,  and  B 
is  defined  by 

(7.2)  eM  -  i  t«.-l  [h^]  .  0  1  i  1. 


Ihe  formula  (7«1)  Is  not  given  In  [33]  but  Is  obtained 
from  the  fomula  given  there  by  analytic  continuation. 

For  finite  values  of  t  we  obtedn  and  Yq 
functions  as  corrections  to  (7<1)*  We  have  [53I 


(7.3)  Xq(h)  -  H(n) 


f (n)  -  u  exp  y'  dt 

■0 


1  -  J  ^  e(t)dt 

0 _ _ 

-I  T  +  2(1  -  J'  ^  e(t)dt^ 


dt 


H(t)([XQ(t)]2 + j 


and 


(7.4)  Yq(m.)  -  H(M,)|g(n)-Hi  exp  J  ~|.M.)dt 


f 

1  ~ 

f  e(t)dt 

0 

T4e( 

c 

^  ^e(t)dt^ 
)  V 

-  1 


dt 


H(t)([XQ(t))2  + 


17? 
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Th*  functions  f  and  g 
(7.5)  t  . 


given  by 


where  P  and  Q  satisfy  the  PredholB  equations 


(7.6) 


P  -  -  L(P)  +  1  +  expl-  T/il  (1  -  nN(-  u)). 


Q  -  L(Q)  +  1  —  exp[-  T/Ul 


1  - 


+  2(1  -f  e(t)dt^ 


with 


(7.7) 


exp 


e(tl 

J 

H(-  1*) - °  l.^U 


dt 


The  linear  Integral  operator  L,  which  we  shall 
not  display «  has  a  continuous  nonnegatlve  kernel.  The 
eqxiatlens  (7.6)  can  be  solved  by  Iteration  to  converge 
xmlformly,  since.  In  this  aom,  we  have  shown  [53]  that 


(7.8) 


IlLlI  ^ 


dt 


The  convergence  Is  very  rapid  except  for  t  near  0. 

A  numerical  program  for  oooqputlng  and  Yq  from 
these  foxnulae  has  a  natural  check.  For  t  ■  0  the 
computations  are  no  different  than  for  t  >  0  and  must 
generate  Xq  and  Yq  to  be  Identically  1.  This  amounts 
to  computing  H  by  (7.1)  and  then  generating  1/R  by 
the  foxmda 
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(7.9)  ir^(n)  -  f  (n)  -  ^  exp 


0 


_ dt _ 

H(t)((Xo(t)]2  +  [^2V/2' 


where  now 


(7.10)  f  =  -  L(f)  +  1  - 


^iwr  ®*p 


I 


0 


e(t) 

%  +  p. 


dt 


We  see  then  that  the  linear  equations  and  linear 
constraints  have  a  unique  solution.  By  the  use  of 
analytic  function  theory  these  are  reduced  to  quadratures 
and  rapidly  convergent  Fredholm  equations. 
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8.  CONCLUSION 

We  have  enirveyed  various  nonlinear  and  linear 
equations  as  possible  substitutes  for  the  Ne\anann  series 
solution  of  the  transport  equation.  We  have  tried  not 
to  forget  the  motivation  for  studying  these  equations « 
that  of  producing  a  slnqpler  conqputatlon  than  the 
Netanann  series. 

The  details  that  we  have  been  able  to  give  are 
possible  because  of  the  existence  of  exact  linear 
equations  among  vdiose  solutions  are  all  the  solutions  to 
the  nonlinear  equations.  We  have  been  able  to  give 
existence  and  iinlqueness  resvilts  as  well  as  results 
relevant  to  the  stability  problem  for  numerical 
conQ)utatlons .  This  wealth  of  detail  Is  primarily  a 
resxilt  of  the  assuo^tlon  of  homogeneity  of  the  atmosphere. 
All  of  the  results  have  been  [52»54l«  or  can  readily  be« 
extended  to  the  anisotropic  but  homogeneous  atmosphere. 

We  will  not  rate  the  various  numerical  methods. 

Since  we  have  made  no  estimates  on  the  rate  of  convergence 
of  Iterative  methods  for  nonlinear  equations^  It  Is  not 
possible  to  coD^are  them  with  the  computations  based  on 
the  linear  Fredholm  equations.  Nor  have  we  given  any 
quantitative  measure  of  the  seriousness  of  the 
Instabilities  discussed  In  Sec.  6.  We  have  shown  that 
Instabilities  put  bounds  on  the  accuracy  of  numerical 
conqputatlons  In  certain  cases.  But  It  may  be  that 
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satlsfactozy  accux^acy  can  be  attained  before  Instability 
sets  In  [  l4] . 

As  we  have  renaxiced  before*  the  linear  singular 
equations  are  a  result  of  the  assumption  of  homogeneity* 
as  clarified  In  Sec.  2.  Numerical  methods  for  solving 
equations  (3*13)  or  equations  (3*1^)  and  (3*15)  perhaps 
have  the  greatest  potential  since  analogous  equations 
can  be  f emulated  for  Inhomogeneous  problems  [lO*19*7^K 
The  detailed  analysis  of  the  linear  equations  for 
homogeneous  problems  and  the  ensuing  results  for  nonlinear 
equations  is  of  value*  if  for  no  other  reason*  in  giving 
indications  of  behavior  to  be  expected  in  a  study  of 
the  more  con^llcated  inhomogeneous  problems. 
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